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Abstract 

Expansion of higher transcendental functions in a small parameter are needed in many ar- 
eas of science. For certain classes of functions this can be achieved by algebraic means. 
These algebraic tools are based on nested sums and can be formulated as algorithms suit- 
able for an implementation on a computer. Examples, such as expansions of generalized 
hypergeometric functions or Appell functions are discussed. As a further application, we 
give the general solution of a two-loop integral, the so-called C-topology, in terms of mul- 
tiple nested sums. In addition, we discuss some important properties of nested sums, in 
particular we show that they satisfy a Hopf algebra. 



1 Introduction 



The expansion of higher transcendental functions [|T], |^] is a common problem occuring in many 
areas of science. It is of particular interest in particle physics in the calculation of higher order 
radiative corrections to scattering amplitudes. There, higher transcendental functions occur 
frequently in formal solutions for specific loop integrals. The necessary expansions of these 
functions are in general a highly non-trivial task. This is particularly true, if the expansions are 
required to a very high order. 

In calculations of higher order radiative corrections classical poly logarithms [|3|], as well as 
Nielsen's generalized polylogarithms [Q] appear. However, this set of functions will not suffice, 
if the number of loops grows, or if several different scales are involved in the problem. Several 
extensions of this class of functions to multiple polylogarithms have been discussed recently [|5|] 
-[§• 

It is the aim of this paper to perform a systematic study of multiple nested sums appearing 
in the expansion of higher transcendental functions around integer values of their indices. To 
that end, we define so called Z-sums, study their algebraic properties and discuss their relation 
to the multiple polylogarithms introduced in the literature [|]] - [||] . We give algorithms to 
solve these multiple nested sums to any order in the expansion parameter e in terms of a given 
basis in Z-sums. All algorithms can be readily implemented on a computer. The Z-sums can be 
considered as certain generalizations of Euler-Zagier sums |T^ or of harmonic sums [|TT]] - 
[ ]T3| ] involving multiple ratios of scales. The latter are known in physics since the calculation of 
higher order Mellin moments of the deep-inelastic structure functions [jlT]], [ [T4l ] - [ [T^ . 

At the same time, our results allow us to investigate higher loop multi-scale integrals occur- 
ing for instance in pertubative corrections to four-particle scattering amplitudes. These integrals 
have received great attention in recent years, mainly motivated by calculations of the next-to- 
next-to-leading order corrections to amplitudes for Bhabha scattering [[T7|], for pp ^ 2 jets [ [TSl ] 
- Jorpp^YY 0] and for light-by-light scattering [^. 

The relevant master integrals at two loops with four external legs have been calculated using 
a variety of techniques. Analytic results were obtained either with the help of Mellin-Barnes 
representations [|23|, |2^] or with differential equations [H, Numerical results were obtained 
by a numerical evaluation of the coefficients of the £-expansion [ p6t ^ . Here, we want to 
advocate a different point of view based on multiple nested sums. As a new result and to 
illustrate our approach, we discuss a specific two-loop integral, the so-called C-topology with 
one leg offshell, which can be reduced for arbitrary powers of the propagators and arbitrary 
dimensions to the aforementioned sums. This is useful for the calculation of the two-loop 
amplitude for e+e^ 3 jets. Some of the techniques presented here have already been used 
in a recent calculation with massive fermions [p8|]. In addition, there exists a wide variety of 
related literature on higher transcendental functions occuring in loop integrals and we can only 
mention a few of them here [ |^ ] - . 

This papers is organized as follows. In the next section we introduce nested sums, show that 
they satisfy an algebra and summarize some important special cases of our definitions. Section|| 
contains the main results of this paper, in particular the algorithms for solving certain classes 
of nested sums. In sec. ^ we give some examples for expansions of generalized hypergeometric 
functions, Appell functions and the Kampe de Feriet function. As an application to higher loop 
multi-scale integrals, we discuss the C-topology in sec. ff^. In Appendix A we show that the 



algebraic structure of nested sums forms a Hopf algebra [ ]36| , pTQ . In Appendix B we briefly 
review the multiple polylogarithms of Goncharov [|^]. 
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2 Definition and properties of nested sums 

We define the Z-sums by 

^ ^ \ 0, n<0, 

n 

Z{n;mi,...,mk;xi,...,Xk) = y-^Z{i-l;m2,...,mk;x2,...,Xk), (1) 
k is called the depth, w = mi + ... +mjt the weight. An equivalent definition is given by 

Z{n;m\,...,mk;x\,...,xk) = V (2) 



In a similar way we define the S-sums by 

S{n) 



n>Ji>!2>...>4>0 



1, n>0, 
0, n<0, 



^(nimi,...,^?^;;^!,...,^^) = Y^^S{i;m2,...,mk;x2,...,Xk). (3) 
Once again an equivalent representation is given by: 



5(n;mi,...,m^;xi,...,Xfe) = £ jW^-'-yk- ^"^^ 

The S-sums are closely related to the Z-sums, the difference being the upper summation bound- 
ary for the nested sums: (« — 1) for Z-sums, i for S-sums. We introduce both Z-sums and 5-sums, 
since some properties are more naturally expressed in terms of Z-sums while others are more 
naturally expressed in terms of S'-sums. We can easily convert from the notation with Z-sums 
to the notation with S-sums and vice versa: 

n il I'l-l '2 

■A- 1 w 1 At 



S{n;mu...;xi,...) = £ t-^ £ Tlj5(j2;m3,...;x3,...) 

,-1=1 n ,2=1 «2 

+S{n;mi+m2,m3,...;xiX2,X3,...), 
Z(n;mi,...;xi,...) = £ £ T-ijZ(/2 - l;m3, ...;x3, ...) 

<l = l '1 12=1 ^2 

-Z{n;mi+m2,m3,...;xiX2,X3,...). (5) 



The first formula allows to convert recursively a S-sum into a Z-sum. The second formula 
yields the conversion from a Z-sum to a S-sum. For example in terms of S-sums the Z-sum 

Z{n;mi,m2,m3,xi,X2,X3) reads 

Z(n;mi,m2,m3,xi,X2,X3) = S'(n;mi,m2,m3,xi,X2,X3) -^(nimi -|-m2,m3,xiX2,X3) 

— 5'(n;mi,m2 + m3,xi,X2X3) +S{n;mi -|-m2 + W3,xiX2X3). 

(6) 
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Figure 1: Sketch of the proof for the multiplication of Z-sums. The sum over the square is 
replaced by the sum over the three regions on the r.h.s. 



Furthermore the Z-sums and the 5-sums obey an algebra. A product of two Z-sums with the 
same upper summation limit can be written in terms of single Z-sums. A straightforward gen- 
eralization of the results given by Vermaseren on the multiplication of harmonic sums yields 

& 

Z{n\mi,...,mk\x\,...,Xk) x Z{n;m[, ...,mi;x[, ...^x'l) 

n 

^Z(/i - l;m2,...,m^;x2,---,-«yt)Z(/i - \;m\,...,m'i,x\,...,x'i) 



n 

+ -j;rZ{i2-Unii,...,mk-,xu...,Xk)Z(i2-Um2,...,m'[;x2,...,Xi) 

'2=1 12 ' 

+ L '^m^w -'^'^2, ■■■,mk;x2, ...,Xk)Z{i - l;m2, ...,m;;4, -,^/)- (7) 

1=1 ' 

Recursive application of eq. (0) leads to single Z-sums. The proof of eq. (0) follows immedi- 
ately from the relation 

n n n i— 1 n 7 — 1 n 

E E = E E + E E + E 

which is sketched in fig. [1|. Note that eq. (0) directly translates into an algorithm for the 
multiplication of two Z-sums. Details on the implementation of this algorithms on a computer 
can be found for example in [[T^. We give an example for the product of two Z-sums: 

Z{n;mi,m2\x\,X2) x Z{n\mj,\x^) = Z(n;mi,m2,m3;jci,JC2,^3) +Z(n;mi,m3,m2;xi,X3,.x;2) 

+Z{n;m'i,m\,m2;x^,x\,X2)+Z{n;m\,m2-\-m^;x\,X2X^) 
+Z{n;mi +m2,m2;xiX2,X2)- (9) 

Note that the product conserves the weight. The Z-sums form actually a Hopf algebra. More 
details can be found in appendix ^ 

The 5-sums also obey an algebra. The basic formula reads 

S{n\mx,...,mk\xx,...,Xk) x S{n\mi,...,mi\Xi,...,Xi) 

E W'^(^'i;"^2, ...,m^;jC2, ...,Xk)S{ii\m\, ...,mi,x\, ...,x'i) 



+ -j;iS{i2\mi,...,mk\xi,...,Xk)S{i2\m'2,...,m!i\x'2,...,x'i) 

12 = 1 ^2 ' 

~ E .^^w '^('"^^2, ...,mk;x2, ...,Xk)S{i;m2, ...,m;;4, (10) 

;'''i-|-m| 
1=1 ' 

Note the minus sign in front of the last term compared to the corresponding formula for Z-sums. 



2.1 Special cases 

Z-sums and 5-sums are generalizations of more known objects. We give here an overview of 
the most important special cases. 

For n = oo the Z-sums are the multiple poly logarithms of Goncharov [Q]: 

Z(oo;mi,...,m^;xi,...,x^) = Li„^,...,„j (x^t, (H) 

For = ... = Xk = I the definition reduces to the Euler-Zagier sums []^, [T0|]: 

Z(n;mi,...,mfc;l,..., 1) = Z„j,...,,„^(n). (12) 

For n = oo and xi = ... = Xk = I the sum is a multiple ^-value 

Z(oo;mi,...,m;t;l,---, 1) = C(mfc,...,mi). (13) 

The 5'-sums reduce for xi = ... = Xk = I (and positive m,) to harmonic sums []T^: 

5(n;mi,...,m;t; I,---, 1) = 5',„i,...,,„^(n). (14) 

The multiple polylogarithms of Goncharov contain as the notation already suggests as subsets 
the classical polylogarithms hinix) [|3|], as well as Nielsen's generalized polylogarithms [Q] 

Sn,p{x) = Lii,...,i,„+i(l^^^,;c), (15) 

p-i 

the harmonic polylogarithms of Remiddi and Vermaseren 

Hmi,...,mif_{p^^ L^m^. , . . . ,mi ( ^1 ; — 7 1^ ; (16) 

k-\ 

and the two-dimensional harmonic polylogarithms introduced recently by Gehrmann and Remiddi 
[§]. The exact connection to the two-dimensional harmonic polylogarithms is shown in ap- 
pendix P together with a brief review of the multiple polylogarithms of Goncharov. Euler- 
Zagier sums and harmonic sums occur in the expansion of Gamma functions: For positive 
integers n we have on the positive side 

r(n+£) =r(i+£)r(R) 

X (l+eZi(n-l)+e2Zii(n-l)+e3Ziii(n-l) + ...+e"-^Zii...i(n-l)). 

(17) 
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Figure 2: The inheritance diagram for Z-sums shows the relations between the various special 
cases. 

On the negative side (again n > 0) we have 



The usefulness of the Z-sums lies in the fact, that they interpolate between Goncharov's multiple 
poly logarithms and Euler-Zagier sums. In addition, the interpolation is compatible with the 
algebra structure. Fig. ^ summarizes the relations between the various special cases. 



In this section we give the detailed algorithms which allow to solve the £-expansion of nested 
transcendental sums in terms of Z-sums or 5-sums defined in eq.(|T]) and eq.(|]), respectively. By 
a transcendental sum we mean a sum over i of finite or infinite summation range involving the 
following objects: 

1. Fractions of the form 



where m is an integer, c a non-negative integer and x a real number. 
2. Ratios of two Gamma functions 



r(-n+l+e) 

r(l+e) (-1)"- 
£ r{n) 



(18) 



3 Algorithms 



r{i + ai+bie) 



(20) 



r{i + a2 + b2ey 
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where ai and a2 are integers. 

3. Z- and 5-sums are also allowed to appear as subsums: 

Z(z + c- l;mi,...;xi,...) or 5(z + c';mi, ...;jci, ...). (21) 
The offsets c and c' are integers. 

4. We further allow these building blocks to occur also with index {n — i), for example as in 

S{n-i;m[,...;x[,...). (22) 
Here n denotes the upper summation limit. 

5. In addition binomials 

"'■ (23) 



/ J i\{n-i)l 
may occur. 

Some examples of sums which can be constructed from these building blocks are 
n y r(/+l+e) - , rji + ae) Tji + be) 

,eiO-+4)3(/+2)8r(/+2 + 3£) r(/+i-c£) ro'+i) • ^ ^ 

There are some simplifications, which can be done immediately: Partial fractioning is used to 
reduce a product for ci ^ C2 

A*-) 1 Aj A^ Aj A^ 



/ + Cl)'"l (Z + C2)'"2 C2-C1 [(Z+Cl)™! (Z + 02)^2-1 (z + Cl )'"! " ^ (z + C2; 



'"2 



(25) 

to terms which involve only the first factor or only the second one, but not both. 

Ratios of two Gamma functions as in eq. ( pO| ) are first reduced to the form r{i + bie) /T{i + b2'£) 
with the help of the identity Y{x + 1) = x T{x). They are then expanded in e using eq. ([T7|). To 
invert the power series which is obtained in the denominator the formula 

(l+£Zi(n-l)+e2Zii(n-l)+£3Ziii(n-l) + ...+£"-iZii...i(n-l))~^ 
= \-zSi{n-\)+£^Sn{n-\)-£^Sn\{n-l) + ... (26) 
is useful to speed up the computation on a computer. 

There are also some basic operations involving Z- or 5-sums. First of all we can easily con- 
vert between the two notations, using eq. (|]). Furthermore we would like to be able to relate 
the Z-sum Z{n + c— 1, ...) to Z{n — 1, ...) or the 5-sum 5(n + c, ...) to S{n, ...), where c > is a 
fixed number. This can easily be done with the help of the following formulae: 

Z(n-fc- \;mi,...\xi,...) 



c-l 

xi—^^ 

5(n + c;mi, ...) 



Z{n-Umi,...;xi,...) + yx\- — -\—Z{n-l+j;m2,...;x2,...), 



S{n;mi,...;xi,...) + yx{ \ S{n + j;m2,...;x2,...). (27) 



1 



The Z- or S-sums appearing in the last term have a reduced depth and the problem can be solved 
recursively. 

Another situation which appears quite often is the product of two sums. If the upper summation 
limits of the two sums differ by some integer c we first synchronize them with the help of eq. 
(p7|). For sums with equal upper summation limit one may use the algebra eq. (0) or eq. ( |T0| ) to 



convert the product into single sums of higher weight. 

Furthermore we can bring Z-sums and S'-sums to a standard form by eliminating letters with 
negative degrees, that is positive powers of i. In general, these cases are easy to handle. We 
illustrate this for 5-sums. We consider S{n; — mi,m2, ...;xi,X2, ...), write out the outermost sum 
of the 5-function and then interchange the order of summation: 



n '2 

A" 



S{n;-mi,m2,...;xi,X2,...) = ^5(/2;m3, ...) z^'^l ■ (28) 



(2 = 1 ^2 ii=(2 

The inner sum can be evaluated for any given weight analytically. Subsequently the outer sum 
can be done with eq. @. If a negative weight occurs inside a sum, eq. ( ^8| ) is applied to the 
subsum starting from the negative weight. 

If a binomial appears in the sum, this sum may be written as a conjugation. To any func- 
tion f{n) of an integer variable n one defines the conjugated function C o f[n) as the following 
sum [0 

Cofin) = -i;(")(-l)Y(/). (29) 

Conjugation satisfies the following two properties 

Col = 1, (30) 
CoCofin) = f{n), (31) 

which can be easily verified. 

We classify four types of transcendental sums, which are dealt with in the algorithms A to 
D. 

1 . Sum over / involving only Z{i — 1 ; . . .) (type A). 

2. Sum over i involving both Z{i — 1; ...) and Z{n — z — 1; ...) (type B). 

3. Sum over i involving S{i; ...) and a binomial (type C). 

4. Sum over i involving both 5(z; . . .) , S{n — z; . . .) and a binomial (type D). 

Many of the algorithms use a recursion. They relate a given problem to a simpler one, either 
with a reduced depth or weight of the Z-sums or 5-sums involved. In these cases we only give 
one step in the recursion. 

The algorithms presented in this paper are all suited for programming in a computer algebra 



system like GiNaC [p8|], FORM [ ]390 or the commercial ones like Mathematica or Maple. Im- 



plementations within the GiNaC framework and in FORM along the lines of ref. [ |T2| ] are in 



preparation or have been published elsewhere pO|]. 
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3.1 Algorithm A 

Here we consider sums of the form 

and show how to reduce them to Z-sums. We assume that all cij and cj are integers, c is assumed 
to be a nonnegative integer and o should be an integer. The upper summation limit is allowed 
to be infinity. 

After expanding the Gamma functions and synchronizing the subsum Z(i + o— l,mi,...) the 
problem is reduced to sums of the form 



x' 



Z(/-l,...) (33) 



with c > 0. It remains to reduce the offset c to zero. If the depth of the subsum is zero, we have 

n yi 1 n i i n 

(i + c)™ ~ x^^ ii + c-\)'"~c^ Jn + cY"' 



The last term contributes only if n is not equal to infinity. If the depth of the subsum is not equal 
to zero, we have 

n i -in i 

y T^^z(i-i,...) = - y 7 — z(/-i,...) 

- E rT-v^42(' - l.'«2, ...) + T^^Z{n - 1, ...). (35) 
Note that the third term only contributes if n is not equal to infinity. Finally we arrive at 

(36) 

which is again a Z-sum. If the upper summation limit n equals infinity this sum yields immedi- 



ately a multiple polylogarithm according to eq. ([11]). In the special case where n equals infinity 



and the subsum is an Euler-Zagier sum we obtain a harmonic polylogarithm according to eq. 

dH). 

3.2 Algorithm B 

Here we consider sums of the form 

x' r{i + ai+bie) r{i + ak + bke) ^ . . 
L 7^~r^F7^~^ rTTV---F7^-^ — — -Z(z + o- l,mi,...,m/,;ci,...,x/) 

y"-' r{n-i + a[+b[e) r{n-i + a'j^,+b'j^,e) 

X 



{n-i + c'Y'' T{n - i + c[ +d[z)"' Y{n - i + c'^, + d[,z) 

xZ{n — i + o' — l,m\, ...,mii,x[, ...,xj/) (37) 
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and show how to reduce them to Z-sums. Here, all Uj, a'j, cj and c'j should be integers, c, c' 
should be nonnegative integers and o, o' should be integers. Note that the upper summation 
limit is (n — l). The upper summation limit should not be infinity. 

Using the expansion of the Gamma functions and the synchronization of the subsums, we im- 
mediately obtain sums of the form 

E— —Z(i— l,mi,...)- —fZin — i— l,m'i, ....). (38) 

Partial fractioning (and a change of the summation index / — > n — / in sums involving the fraction 
with {n — i + c')) reduces these sums further to sums of the type 



Z{i — l,mi, ...)Z{n — / — l,mj, ....). (39) 



If the depth of Z{n — / — l,mj, ....) is zero, we have a sum of type A with upper summation 
index in — l): 



Z(/-l,mi,...). (40) 



Otherwise we can rewrite eq. ( |39D as 

n-l 



E 



7-1 yi yi-' 

E— —Zii— l,mi,...) — ^- rZij — i— Lmo,....! 
(i + c)™ ^ ' ' 2. ; 



(41) 



and use recursion. The inner sum is again of type B, but with a reduced depth, such that the 
recursion will finally terminate. 

3.3 Algorithm C 

Here we consider sums of the form 

n\,_^.i ^' T{i + ai+biz) T{i-\-ak-\-bkZ) 

i=l 

xS{i + o,m]_,...,mi,Xi,...,xi), (42) 

where aj and cj are integers, c is a nonnegative integer and o is an integer. The upper summation 
limit should not be infinity. These sums cannot be reduced to Z-sums with upper summation 
limit n alone. However, they can be reduced to Z-sums with upper summation limit n and mul- 
tiple poly logarithms (which are Z-sums to infinity). 

Again, we expand the Gamma functions and synchronize the subsum. It is therefore sufficient 
to consider sums of the form 

t(:)(-'Yj^S(,,...) (43) 



(=1 
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with c > 0. To reduce the offset c to zero we rewrite the sum as 

.i) )(_,). ,...). (44) 

Repeated application of the above relation yields sums of the form 

-i{")(-nfs{i,-). (45) 
If m is negative we rewrite eq. as 



We can therefore assume that m is a non-negative number. Furthermore, due to eq. ( PS] ) we may 
assume that in the 5-sum 5(z;mi, ...) no mj is negative and that if some mj is zero, then 
the corresponding xj is not equal to 1. The sum 5(z, ...) is then rewritten as 

S{i;mi,...,mk;xi,...,Xk) = S{N;mi,...,mk\xi,...,Xk) 
-S{N;m2,...,mk\X2,...,Xk)-x\ ^ 

\h=i+l 



N N 

+5(A^;m3,...,m^;.JC3,...,Jc^) X I 



z ' z'"2 

!2=(i + l 1 2 

A' A' N i\ il ik 

Al Ao A;. 



-••• + (-l)M I I - I ^%...^A. (47) 
The proof of eq. ( P7| ) is not too complicated and consists in repeated application of the identity 

n i N i N N N N 

(=1 7=1 ;=1 7=1 7=)i+l 7=1 7:=(i+l ;■=;■+ 1 

Eq. holds for any and in particular we may take = 0° in the end. Each term is then a 
product of a 5-sum at infinity and a sum of a new type. The 5-sum at infinity is converted to a 
Z-sum at infinity and expressed in terms of multiple poly logarithms. We now deal with sums of 
the form 

" f n \ ■ x|) ^ ^ ^ x'^ X2 x^l 

~E( E E-. E (49) 

We introduce raising and lowering operators as follows: 

(x+)"M = i-ln-W, 

X 

dx 



x--/(x) = xj-J{x). (50) 
11 



It is understood that in the second line only functions which are integrable atx = are consid- 
ered. We see that x~ is the inverse to x"*", e.g. x~x"'" = id. However x+x^ = id holds only if 
applied to non-trivial sums. For the trivial sum we have x^x^Z{n) = 0. 
With the help of the raising operators eq. (^9|) may be rewritten as 



(xk+)™^ (xk-i+)™^-' ... (X1+)"" {xo+r (-1) I ( ■ ) (-^0 



N N N 

X I I - I (51) 

It may happen that some x, 's are equal to one. In this case we first calculate the sum for arbitrary 
jc/'s and take then the limit Xi ^ 1. Some care has to be taken for the double limit x ^ I and 
N oo. The order is as follows: First all limits x — > 1 are taken, then the limit N in eq. 
( p7| ) is performed. 



The sums in eq. ( pTD can be performed with the help of the geometric series 

^ XX 

y ^' = (52) 

It is evident that if we don't have to take the limit x ^ 1 we can immediately neglect the second 
term. Also in the case x = I the second term can be neglected. It gives rise to terms of the form 

(-Tj^,^ = Lyr (53) 



On the r.h.s the limit x ^ 1 may safely be performed and the resulting sum gives a vanishing 
contribution in the limit N ^ 



Performing the sums in eq. (|T]) we therefore only have to consider expressions of the form 

l-Xk\-Xk-lXk l-Xl...Xk 

x[l-{l-xoXi...Xkr]. (54) 

We then perform succesivly the integrations corresponding to the raising operators. The basic 
formulae are: 

x+[i-(i-W] = £-[i- (1-^1^2; 



x!T^^[i-(i-w2r] = -(i-^ori-f-r^) 

1 -X1X2 I \l-xoJ 



n 1 / 1 \ ' 

1 - (1 -XoXlX2y 



+(i-(i-^or)i;^ 



(=1 



+x+--^^(^ix2f (i-(i-xor), 
1 —X1X2 

1 —X1X2 n ^ I 



+ (x+)-p^^(xiX2f . (55) 
^ ' I -X1X2 
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We use the first fonnula to reduce mq in eq. (Q) to zero: 

n 1 

El / / -\-\mi / -^\mo-l Xk Xjc—lXk Xl...Xk 



1 - (1 -X0Xi...Xk) 



(56) 



In the following we may therefore assume mo = in eq. (|5^. If mq = 0, mi > and xqj^ I we 
obtain for eq. 



-(.-.„)«eKt^ 



Xo 



l-Xkl-Xk-\Xk" l-X2...Xk 



l-{\-{xQX\)...XkY 



N J N 



A' yh 'k 

At A;, 



,■=1 ' i2=i+l ik=ik-l + i 2 'fc 



(57) 



Here we neglected terms of the form eq. (|53|). In the case mo = 0, mi > and jco = 1 we use 
the third formula of eq. (|55|). Again we may neglect contributions of the form eq. (|53|). Doing 
so we obtain 



1 



' ^ ^ ' 1 -Xk 1 



mi-l .X:^ Xk-\Xk X2---Xk 



Xkl-Xk-lXk l-X2...Xk 
N jl Jk 



[1 -(1 -XI. ..Xk) 



^2 



N N 



(58) 



It remains to treat the last term to complete the recursion. The last term introduces a sum of the 
type 



N N Jl Jk 

E ■ E 



il=n+l ik=ik-l + i 1 

Using the inverse formula to eq. (^7|) 



^1 

■mi ••• -ink ■ 
'1 



(59) 



N 



N il Jk 

-^k 



^2 ■ ■ ■ ^2 ■ ■ ■ -"^/t 

ii=n+l ik=ik-l + '^ 'l '■k 

= {-lfS{n\mi,...mk\xx,...,Xk) - {-\fS{N\mx,...mk\xx,...,Xk) 
+{~lfS{N\m2,...mk;x2,...,Xk) 

+(-l)^S(Mm,;xfc) - E 

zi=;3+l ik-\=ik-2+'^ 



(60) 



this sum is easily related to 5-sums with upper summation limit n and (after taking the limit 
N ^ °o)io multiple poly logarithms. 
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If mo = 0, mi = and xi 7^ 1 we rewrite eq. ( p4| ) as 
1 



Xk X3...Xk {xiX2)x3...Xk 

•••I 3; I 21 j l-Xk"l-X3...Xkl-{xiX2)x3...Xk 



X 



1 - (1 -Xo{xiX2)x3...Xky 

X\ I +\int I +\in2 Xk 



+ Y^W)'"'...(4)"'y- 



X2X3...Xk 



■Xk \-X2X3...Xk 



1 - (1 - {xQXi)x2X3...Xky 



(61) 



The case mo = 0, mi = and x\ = \ has to be excluded. However, with an appropriate choice 



of the standard form for S'-sums (c.f. eq. ^8|) this case never occurs 



3.4 Algorithm D 

Here we consider sums of the form 

f^\i J {i + c)"'T{i + Ci+dizy"T{i + Ck + dkt) 

x5'(/ + o,mi, ...,m/,.x;i, 

T{n-i + a[+b[z) Y{n-i + a',^, + b[,t) 



X 



(n - / + c')""' T{n - i + c[ +d[z)"' Y{n - i + c[, + d[,z) 
xS{n — i + o\m[, ...,m\t,x\, ...,x'ii). 



(62) 



Here, all a ,, a'-, c, and c' are integers, c, c' , are nonnegative integers and o, o' are integers. Note 
that the upper summation limit is (n — 1). The upper summation limit should not be infinity. 
As in the case of sums of type C, we cannot relate these sums to Z-sums with upper summation 
limit {n — l) alone, but we can reduce them to Z-sums with upper summation limit (« — 1) and 
multiple poly logarithms (which are Z-sums to infinity). 

After the expansion of the Gamma functions and the synchronization of the subsums we have 
sums of the form 



1=1 



(-1)'77 



(/ + c)' 



-5(?,mi,...)- 



n — i + C 



—^S{n-i,m[, 



(63) 



Partial fractioning leads to 

-E 

i=l 



X 



(64) 



with c > 0. In order to reduce the offset c to zero one rewrites eq. ( |64{ ) as 



1 



1 



X J n + 1 



-V 



i=l 



n+1 
i 



X' 



{i + c-iy 



-iS{i — 1, ...)5(n + 1 — z,mi, 



We arrive at sums of the form 



-Si" 



(65) 



(66) 



14 



If the depth of S{n — z, m'^, ....) is zero, we have a sum of type C: 



n-l 



(=1 



L[ i j(-l)'^S(/,mi,...) + ^S(n,mi,...). 



(67) 



Otherwise, we first reduce m to zero. For m > we rewrite eq. (g6|) as 



n 1 



7=1 



J 



7-1 



n 1 



(68) 



For m < we rewrite eq. (^51) as 



n-2 
;=1 



n-l 
i 



+n.x:5(l,mi, ...)5(n — \,m\, ...). 
Having reduced m to zero we arrive at sums of the form 



(69) 



(70) 



For jc 7^ 1 we obtain after some algebra 



^ 1 / 1 N-' 



+ 



;=1 



(-1)1 • (-i)'^'^O' 



(n-iri ^ 



(71) 



where the original problem is reduced to one of the same type but with lower weight. In the 
case x = I the r.h.s of eq. ([71]) reduces to 



1 



-(-1) I ( ■ J (-l)'^5(/,m2,...)5(n-z,m;,....) 



/=l 



1 c/ • / 

n — ipi ^ 



(72) 



Again, the original problem is reduced to one of the same type but with lower weight. The 
above algorithm yields thus a recursion to treat sums of type D. 
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4 Applications 



The algorithms given in this paper can be used for the expansion of higher transcendental 
functions around integer values of their indices, where the expansion parameter occurs in the 
Pochhammer symbols. In this section we give a few examples. Additionally, we illustrate the 
applicability of the algorithms for nested sums to the calculation of loop integrals, in particular 
to integrals with several scales. As an example, we discuss the C-topology and show that the 
integral can be written as a nested sum of the type previously discussed. 

4.1 Generalized hypergeometric functions 

The generalized hypergeometric functions are defined by [[l|, 

j+iFj{au...,aj+i;bi,...,bj;x) = 2^ ^j—. jj—. r, (73) 

where = Y{n + a) /T{a) denotes a Pochhammer symbol. These functions can be rewritten 
as 

j^iFj{a\,...,aj+\;b\,...,bj\x) 

^ Y{bi)...Y{bj) - TO' + ai) T{i + aj) T{i + aj^,) 

r(ai)...r(ay+i),ti ni+bi)-T{i+bj) r(/+i) 

and fall therefore into the category of transcendental sums of type A. We give a few examples 
obtained using the algorithms given in sec. p7T| : 

2Fi(a8,te;l -cz;x) = l+abhi2{x)£^ + ab{c'Li3{x) + {a + b + c) Si,2W)£^ + 0(£'^), 
2Fi(l, -£; 1 -£;x) = l+ln(l-x)£- Li2(;c)£^- Li3(jc)£^- Li4(;c)£'^- Li5(x)£^ 

-Li6(;c)£^- Li7(jc)£^ + 0(£^), (75) 

3F2(-2£, -2£, 1 - £; 1 - 2£, 1 - 2£;jc) = 1 + 4 Li2(x)£2 + (12 Li3(jc) - 45i,2(^)) £■' 

+ (32 U4{x) +45i,3(x) - 1 252,2 (■^) ) £^ + (80 Li5(jc) - 45i,4(x) + 1 2^2,3 (x) - 32^3,2 (x) ) £^ 
+0(£^), (76) 

which all agree with known results in the literature ^ . 

4.2 Appell functions 

The first Appell function is defined by |^] 

Fi{a,biM\c;xuX2) = L 2^ ( \ i f ^^^^ 

mi=0m2=0 \C)my+m2 m\.m2. 

It can be rewritten as 

Fi{a,b\M\c\x\,X2) 

r(c) - , r(/+a)r(/+z.i) r(c) - , r(/+a)r(/+ZP2) 
r(a)r(Z7i) ' r(/+c)r(/+i) ^T{a)T{b2)ti ' r(/+c)r(/+i) 
r(c) - r(n+a) "^i , r(/+Z7i) „_,. r(n-/+z.2) 



'T{a)T{bi)T{b2)ti^{n + c)f^, ' T{i+\) ^ T{n-i+\)- 
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The inner sum of the last term is of type B. The first Appell function can therefore be expanded 
with the help of the algorithms A and B. 



The second Appell function is defined by [ p3| , ^ 



F2{a,bi,b2;ci,C2;xi,X2) 



oo oo 

E E 



(^)rai+m2(^l)mi (^2)m2 '''•1 -^2 
Om^O (Cl)mi(c2)m2 mi!m2!' 



mi- 



(79) 



It can be rewritten as 



F2{a,bub2;ci,C2;xuX2) 



r(ci) y ^ r{i+a)rii+bi) r(c2) y ,r(^+a)r(H-^ 
^ r(a)r(Z7i) ^tt""' r(/+ci)r(/+ 1) ^ r(a)r(Z72) ,tt'''r(/+c2)r(/+ 1) 

r(ci)r(c2) ^ r{n+a) 
~r(a)r(z>i)r(z>2)„tir(n + i) 



(=1 



(-1)1 ■ (-l)'(-^l) 



r(/+Z7i) ^^„_, r(n-/+Z72) 
r(/+ci)'^2 r(n-/+c2) 



(80) 



The inner sum of the last term is of type D. The second Appell function can therefore be ex- 
panded with the help of the algorithms A to D. As an example we give 

^2(1, !,£;!+£, 1 -£;a:,3;) = + (21n(l -x) -ln(l -x-y))e 



1 —X 1 —X 



1 



1-JC 

2Lii,i 
1 



2U2{x) + 2U2{y)-U2{x+y)+4So_2{x)+So,2{x + y)+Ui^i (^,x+y] 

\x + y J 



x + y 



,x + y -2Lii,i 



x + y 



-,x 



Li3 (x + y) -2 Us (x) + 5o,3 {x+y)- 85o,3 (x) - Si ,2 (x + y) - 45i ,2 (x) + 45i ,2 (y) 

y 



l-x 

Hi,2{x+y) -4Hi^2{x) +2Ui^2 
x + y 



+2 Li, 2 



X +2Li2,i 



-4Li2,ig,x)+4Lii,i,i(^ 



x + y 
X x + y 



+4Li> 1, 1 



+ Lii , 1 



x + y 



x + y X 



x + y +2Lii,i,i 1 



-—,x + y - Lii,2 
x + y ) 

x + y^ +3 Li2,i 
,x\ +4Lii 1,1 



x + y 



.x+y 



y 



x + y 
x + y 



y 



x + y 



x + y -2Lii,i, 



x + yj +2 Li2,i 
hx^ 

y x+y 



x + y 



X 



-,x 



x + y X 



,x 



x + y 



l,x + y -2Lii,i,i 



-2Lii.ij ( ,— — ,jc + 3; 

X x + y 



x + y 
e^ + 0{e^) 



A, x + y -2Lii,i,i 1, 



x + y 



(81) 



After taking into account a typo in eq. (A.47) of ref. [ |44| ] this result agrees up to order £ with 
the one obtained along the lines of ref. [ p4| ] . Multiple poly logarithms of low weight can be 
expressed as products of classical polylogarithms and the result of the expansion in eq. ( |8T1 ) 
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can be simplified Q. However, we present here the output of our algorithm. Our algorithm 
gives the result as a linear combination of polylogarithms and does not generate products of 
poly logarithms. 



4.3 Kampe de Feriet function 

The Kampe de Feriet function is defined by 

Si[aua2,bi;c,ci;xuX2) =1^1^ -p. ; — n j- — j-- (o2) 

It can be rewritten as 

Si{a\,a2,b\;c,ci\x\,X2) 

r(c)r(ci) y , r(/+ai)ro-+a2)r(/+z>i) 
r(ai)r(a2)r(z>i),ti'^' r(/+c)r(/+ci)r(/+i) 
r(c) y i r(/+ai)r(/+a2) 
+r(ai)r(a2),tt''' r(/+c)r(/+i) 

r(c)r(ci) y „ r{n + ai)r{n + a2) 
r(ai)r(a2)r(Z>i),^/2 r(n + c)r(n + l) 

^(-.)l(:)(-.r(-|)'R^. 

The inner sum of the last term is of type C. The Kampe de Feriet function can therefore be 
expanded with the help of the algorithms A and C. 



4.4 The C-topology 

Here we study the C-topology with one massive external leg and arbitrary powers of the prop- 
agators and arbitray dimensions, which can be solved using the algorithms given in this paper. 
The importance of this result lies in the fact that one avoids having to solve a system of equations 
obtained from partial integration or Lorentz invariance. Solving such a system becomes quite 
difficult if one external leg is massive. The result obtained here is thus a useful ingredient for the 
calculation of the two-loop amplitudes with one massive external leg, such as e+e^ 3 jets. 
The second C-topology, where the massive leg is attached to the other comer, as well as all 
simpler topologies, can be obtained along the same lines. The two-loop C-topology with one 
massive external leg is defined by 



d^ki r d^h 11111 

h^J m«/2 (-^2)Vi ^_^2)V2 ^_^2)V3 ^_^2)V4 ^_^2)V5 



(84) 



with 



h = h-P2, 

k4 = h-pn-i- (85) 



' We would like to thank the referee for pointing this out. 
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Figure 3: The C-topology reduces to a triangle with three external masses and an additional 
integration over the Feynman parameter a. 



Fig. H shows the corresponding Feynman diagram. We first perform the /s -integration. Com- 
bining with Feynman parameters first /| and /| and then the resulting propagator with Z| we 
obtain: 



r(V235 - m + e) r( -V23 + m - e)r( -V5 + m - e) /■ , V2-1-V3-1 
:^(V2)r(V3)r(V5) r(-V235 +2m-2£) J 

(86) 



r(V2)r(V3)r(V5) r(-V235 +2m-2£) 

u 

d^kl 1 1 1 



As a short hand notation we used D = Im — 2z,d=\—a and V235 = V2 + V3 + V5 . The second 
line is a one-loop triangle with three external masses. The solution for this one-loop integral 
with arbitrary powers of the propagators and arbitrary dimensions is known [|45|]. We use the 
solution given in [M and perform the remaining integration. We obtain 



r(2m - 2e - Vi235)r( 1 + V1235 - 2m + 2e)r(2m - 2e - V2345)r( 1 + V2345 - 2m + 2e) 

r(vi)r(v2)r(v3)r(v4)r(v5)r(3m - be - V12345) 

^^ r(m-e-V5)r(m-e-V23) , ^2m-2e-Vl2345 y 

r(2m-2£-V235) ^ '"'^ hh''-'^- 
r(/i + V3)r(/2 + V2)r(/i + /2 - 2m + 2£ + Vi2345)r(/i + /2 - m + £ + V235) 



X 
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r(/i + 1 - 2m + 2£ + Vi235)r(/2 + 1 - 2m + 2£ + V2345)r(/i + ii + V23) 

2m-2e-Vi235 

r(;i + 2m - 2£ - Vi25)r(/2 + V2)r(/i + i2 + V4)r(/i + /2 + m - £ - Vi) 
r(/i -f 1 +2m-2£-Vi235)r(/2 + l - 2m-F2£-FV2345)r(/i -f/2 + 2m-2£-v 

,2m— 2£— V2345 

r(/i + V3)r(/2 + 2m - 2£ - V345)r(/i -I- i2 + Vl)r(/i /2 + m - £ - V4) 
^ r(?i + 1 - 2m + 2£ + Vi235)r(/2 + 1 + 2m - 2£ - V2345)r(/i + ii + 2m - 2£ - V45) 

2m-2£-Vi235 2m-2£-V2345 r(/i + 2m - 2£ - Vi25 )r(/2 + 2m - 2£ - V345) 



X 



X2 



T{h + 1 + 2m - 2£ - Vi235)r(/2 + 1 + 2m - 2£ - V2345: 
r(/i + 12 + 2m-2e- V235)r(/i + 12 + 3m - 3£ - v 12345) 



r(/i + Z2 + 4m - 4£ - V 12345 - V5) 



(87) 



where we setxi = (— 5'i2)/(— ■^123) andx2 = (— 5'23)/(— ■^123)- Changing the summation indices 
as in the case of the second Appell function yields a sum of type D. For specific (integer) values 
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of V, and m this expression can be expanded in £ with the algorithm given in sec. ^. This is 
a new and useful result. Up to now this integral has only been known for m = 2 and the sets 
(1, 1, 1, 1, 1) and (1, 1, 1, 1,2) for (vi, V2,V3, V4, V5). We have verified that our result agrees in 
these two specific cases with the ones given by Gehrmann and Remiddi [|]]. 



5 Conclusions 

In this paper we studied some algebraic properties of nested sums. Based on these properties we 
developed a number of algorithms which can be used to expand a certain class of mathematical 
functions. All presented algorithms are suitable for the implementation on a computer. These 
algorithms allow the evaluation of integrals occuring in high-energy physics. As an application 
we have shown how the two-loop C-topology can be evaluated for arbitrary powers of the prop- 
agators and arbitrary dimensions. Furthermore we have shown that the nested sums satisfy a 
Hopf algebra and established the connection with the Hopf algebra of Kreimer. 
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A The Hopf algebra of Z-sums 



In this section we show that the Z-sums form a Hopf algebra [ pq , pTQ . It is sufficient to demon- 
strate that the Z-sums form a quasi-shuffie algebra. A general theorem guarantees then, 
that they also form a Hopf algebra. We also discuss the connection with the Hopf algebra of 
Kreimer [BTl]. 



Before we start, we have to introduce some notation. We call a pair [mfXj) a letter and the 
set of all letters the alphabet A. We further call mj the degree of the letter (nij^Xj). On the 
alphabet A we define a multiplication 

{mi,xi) ■ {m2,X2) = (mi +m2,xiX2), (88) 

e.g. the x/s are multiplied and the degrees are added. As a short-hand notation we will in the 
following denote a letter just by Xj = {mfXj) . A word is an ordered sequence of letters, e.g. 

W = Xi,X2,...,X^ (89) 

We denote the word of length zero by e. The sums defined in ([T]) are therefore completely 
specified by the upper summation limit n and a word W . In particular for any positive n the sum 
corresponding to the empty word e equals 1. A quasi-shuffle algebra SI on the vectorspace of 
words is defined by [ p6| ] 

eoW = Woe = W, 

{XuWx)o{X2,W2) = Xu{Wio{X2,W2))+X2,{{XuWi)oW2) 

+{X^-X2),{WioW2). (90) 
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Note that "•" denotes multiplication of letters as defined in eq. (^8|), whereas "o" denotes the 
product in the algebra Jl, recursively defined in eq. (^D|). We observe that the formula for the 
multiplication of Z-sums eq. (0) is identical to eq. (|90|). The Z-sums form therefore a quasi- 
shuffle algebra. 



We now discuss the connection with the Hopf algebra of Kreimer [|47|]. Kreimer showed that 
the process of renormalization of UV-divergences occuring in quantum field theories can be 
formulated in terms of a Hopf algebra structure. We first recall the properties of an algebra and 
a coalgebra: An algebra has a unit and a multiplication, whereas a coalgebra has a counit and a 
comultiplication. A Hopf algebra is an algebra and a coalgebra at the same time, such that the 
two structures are compatible with each other. In addition there is an antipode. We show that 
the coalgebra structure of Z-sums is identical to the coalgebra structure of the Hopf algebra of 
Kreimer. To this aim we introduce the explicit definitions of the counit, the coproduct and the 
antipode. It is convenient to phrase the coalgebra structure in terms of rooted trees. A Z-sums 
can be represented as rooted trees without any sidebranchings. As a concrete example we write 
down the pictorial representation of a sum of depth three : 



Z(n;mi,m2,m3;xi,X2,X3] 



/l = li2=l;3 = l 1 ■ 



Xl • 

X2» 
Xsi 



(91) 



The pictorial representation views a Z-sum as a rooted tree without any sidebranchings. The 
outermost sum corresponds to the root. By convention, the root is always drawn on the top. 

Trees with sidebranchings are given by nested sums with more than one subsum, for exam- 
ple: 




Of course, due to the multiplication formula, trees with sidebranchings can always be reduced 
to trees without any sidebranchings: 

" x' 

Y^^-^^i.^-'^'^"^2,X2)Z{i-l\m3\X3) = 

Z(n;mi,m2,m3;xi,X2,X3) +Z(/7;mi,m3,m2;xi,X3,X2) 

+Z{n;mi,m2 + mT,\xi,X2X-i). (93) 

The coalgebra structure is formulated in terms of rooted trees (e.g. there is no need to convert 
rooted trees to a basis of rooted trees without sidebranchings). We first introduce some notation 
how to manipulate rooted trees. We adopt the notation of Kreimer and Connes An 
elementary cut of a rooted tree is a cut at a single chosen edge. An admissible cut is any 
assignment of elementary cuts to a rooted tree such that any path from any vertex of the tree to 
the root has at most one elementary cut. An admissible cut maps a tree ? to a monomial in trees 
tio ...otk^i. Note that precisely one of these subtrees tj will contain the root of t. We denote 
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this distinguished tree by R^{t), and the monomial delivered by the k other factors by P^{t). 
The counit e is given by 

e{e) = 1, 

e{t) = 0, tj^e. (94) 
The coproduct A is defined by the equations 

A(e) = 

A{t) = e(^t + t^e+ P'^{t)(^R'^{t), 

adm. cuts C of / 

A{ti0...otk) = A(fi)(o®o)...(o®o)A(ffc). (95) 
The antipode S is given by 

5(e) = e, 

S{t) = -t- £ s{p^{t))oRC(t), 

adm. cuts C of I 

S{ho...otk) = Sih)o...os{tk). (96) 
The proof that these definitions yield a Hopf algebra can be found in [BH] and is not repeated 



here. The Hopf algebra of Kreimer and Connes [ pTj p8| ], which emerged in the context of 
renormalization of UV-divergences, has a slightly different algebra structure. There the algebra 
is generated by rooted trees. In this algebra a product of two rooted trees is not necessarily a 
rooted tree again. However, the coalgebra structures are identical, which is a remarkable obser- 
vation. 

Let us give some examples for the coproduct and the antipode for Z-sums: 

AZ(«;mi;jci) = e®Z(n;mi;xi) +Z(n;mi;xi) 
AZ(n;mi,m2;JCi,JC2) = e®Z(n;mi,m2;xi,X2) +Z(n;mi,m2;JCi,JC2) ® e 

+Z{n;m2;x2) ®Z{n;mi;xi), (97) 

SZ{n;mi;xi) = —Z{n;mi;xi), 
SZ{n\mi,m2\x\,X2) = Z{n\m2,mi\X2-,x\) +Z{n\mi+m2\x\X2). (98) 

B Review of Goncharov's multiple poly logarithms 

At the end of the day we express our results in terms of Goncharov's multiple poly logarithms. 
They form therefore an important specialization of nested sums and we review therefore some 
of their properties. After the introduction by Gonacharov [Q| they have been extensively studied 
by Borwein et al. [§]. They use a different notation which is related to the one of Goncharov by 

Li„„...,„,(x,,...,xi)=^f bj = ^ . (99) 



b\,...M y' xiX2...Xj 
Most of the material reviewed in this section is based on the work of Borwein et al. 
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B.l Integral representations 

We first define the notation for iterated integrals 

dt dt f dtn f dtn-i f dti 
0...0 = / / X ... X / . (100) 

Qn—t a\—t J a„—tfiJ dn-l^^n-l J ^^1 ~ h 



We further use the following short hand notation: 

A A 

I dt \ dt f dt dt dt ,,r.,^ 

— o = / — o... — o . (101) 

\t J a — t J J ^ t ^ a — t 

m times 

The integral representation for \lmi„...,mi {^k^---:^\) reads: 



X2X2,...Xk — t 



o — o o ... o — o 



t J X2,...Xk — t \t J 

In the notation of Borwein et al. this representation reads 

dt dt 



(102) 



Lim^,...,mi(^yt,...,^l) = (-1)^ / ( — o 



t I t-b] 



dt dt fdt \"''~^ dt 



o — o 



t I t — b2 \ t J t — bk 



where the b/s are related to the x/s as in eq. (P9|). Changing the integration variables according 
to t ^ 1 —t yields the dual integral representation: 



1 ■ • mt-l 



Li„„...,„,(x,,...,xi) = (-i)Vr:^°(T^°) 





dt ( dt \ dt ( dt ^'"'"^ 



O 0...0 o 



i-t>k-i-t \\-t J \-bi-t \\-t 



(104) 



In addition to these weight-dimensional integral representations there is also a depth-dimensional 
integral representation: 



1 [ dti (ln?i 



1 (mi).. A (nik) J 



\tn] —1 



r{mi)...r{mk)J ti ^-1 

.,^^(!5f^....></>(l^^|r^^ ao5, 
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B.2 The shuffle algebra 

Instead of specifying a multiple polylogarithm by the x/s and m/s, we may denote the function 
by a single string 

(ai,a2,...,aw) = {0,...,0,bu0,...,0,b2,...,0,...,0,bk) , (106) 
where {nij — 1) zeros preceed bj. Defining further 

a(ai) = (107) 

t-Ui 

allows us to rewrite the integral representation eq. ( |103D as 

1 

Li„,,...,„i(xfc,...,xi) = (-l)^y"a(ai)o...of2(a^). (108) 



From the iterated integral representation one deduces a second algebra structure with multipli- 
cation given by 

^hni^,...,mi {^ki ■■■■i^l) ^ Lim^._|_;,...,nn._|_i {^k+h ■■■^■^k+l) 
1 1 



1 

= (-1)'+' E /f2(a^(i))o...oa(a^(,,^+,,)), (109) 

shuffle-'^ 

where Wk = m\-\- ... +mk, wi = m^t+i + ••• +i^k+i and the sum is over all permutations, which 
preserve the relative order of the strings t2(ai)...t2(awj and Q.{awi,+i) ■ ■ -^{cLwh+wi) ■ 

B.3 Reduction to simpler functions 

The multiple polylogarithms contain a variety of other functions as a subset. We start with depth 
one. As the notation already suggests, the multiple polylogarithms are in this case identical to 
the classical polylogarithms, e.g. 

Lio(jc) = ^^, Lii(x) = -ln(l-x) (110) 
1 —X 

and 

X 

UJx)= fdt^^^. (Ill) 







Nielsen's generalized polylogarithms [^, defined through 
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are related to the multiple polylogarithms by 

Sn,pix) = LiL...j,„+i(l^^^,x), (113) 

p-l 

where (p — I) one's occur before n + l and x. The harmonic polylogarithms of Remiddi and 
Vermaseren [0] are related to the multiple polylogarithms for positive indices as 

Hmi,...,mi{p^^ ^hni^,...,m\ Q ; —7 ^ 7-^) ■ (114) 

k-\ 

The harmonic polylogarithms are defined recursively through 

//o(x) =ln(jc), Hi{x) = -\n{\-x), //_i(jc) = ln(l +jc) , (115) 

and 

«... - 



X 

H±l^m2,...,mk = J dtf±i{t)H,n2,...,mki^), (116) 


where the fractions fo{x), fi {x) and f-i{x) are given by 

/oW = -, /i(^) = -^, f_,{x) = -^. (117) 

Recently Gehrmann and Remiddi [||] extended the harmonic polylogarithms to two-dimensional 
harmonic polylogarithms (2dHPL) by extending the fractions to 

/(z,x) = ^, /(l-z,x) = -J . (118) 

z+x l—z—x 

From the integral representation eq. ( |lU2p it is clear that the 2dHPL are a subset of Goncharov's 
multiple polylogarithms. If we identitfy x = x\X2 x ... x ;cyt we have 

/d-tf{z,t)lAffi, [ Xj(^, ...,X2f 1 = ~Limt,...,mi,l \Xjc,....,X2, , ) , 

\ X2...XkJ \ X2...Xk Z 



/dtf\\ Z, ?)Limi.,....mi \ X]i^ ■■■^X2-, 1 = ^^m],,...,m^A \ X.]i^i ■■■^X,2^ , -; 1 • 

V X2...XkJ V X2...Xk l-Zj 





X 



(119) 
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